Excitonic and vibronic spectra of Frenkel excitons in a two-dimensional simple lattice 
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Excitonic and vibronic spectra of Frenkel excitons (FEs) in a two-dimensional (2D) lattice with 
one molecule per unit cell have been studied and their manifestation in the linear absorption is 
simulated. We use the Green function formalism, the vibronic approach (see Lalov and Zhelyazkov 
[Phys. Rev. B 75, 245435 (2007)]), and the nearest-neighbor approximation to find expressions of 
the linear absorption lineshape in closed form (in terms of the elliptic integrals) for the following 2D 
models: (a) vibronic spectra of polyacenes (naphthalene, anthracene, tetracene); (b) vibronic spectra 
of a simple hexagonal lattice. The two 2D models include both linear and quadratic FE-phonon 
coupling. Our simulations concern the excitonic density of state (DOS), and also the position and 
lineshape of vibronic spectra (FE plus one phonon, FE plus two phonons). The positions of many- 
particle (MP-unbound) FE-phonon states, as well as the impact of the Van Hove singularities on 
the linear absorption have been established by using typical values of the excitonic and vibrational 
parameters. In the case of a simple hexagonal lattice the following types of FEs have been considered: 
(i) non-degenerate FEs whose transition dipole moment is perpendicular to the plane of the lattice, 
and (ii) degenerate FEs with transition dipole moments parallel to the layer. We found a cumulative 
impact of the linear and quadratic FE-phonon coupling on the positions of vibronic maxima in the 
case (ii), and a compensating impact in the case (i). 

PACS numbers: 71.20.-b, 71.35.-y, 71.35. Aa, 71.35. Cc 



I. INTRODUCTION 

Frenkel excitons (FEs) and their vibronics have been 
studied for several decades. ^'^ The theoretical basis of 
the vibronic studies and their connection with spectro- 
scopic data have been established in papers,'^"'' as well 
as in a couple of reviews and books. ^"^^ The studies of 
charge transfer systems enlarge the applicability of the 
exciton theory towards the concept of charge transfer ex- 
citons (CTEs) and their vibronics. The coupling be- 
tween FE, CTEs, and phonons complicates the vibronic 
spectra manifesting themselves in various linear and non- 
linear phenomena. Several mechanisms of this coupling 
combined with the intermolecular transfer of quasipar- 
ticles create bound (one-particle) exciton-phonon states 
and unbound many-particle (MP) states with rather dif- 
ferent lineshapes. The models of vibronic spectra need 
productive methods of calculations and interpretation of 
the spectral pictures. 

In this paper, we calculate the vibronic spectra of a 
two-dimensional (2D) lattice with one molecule per unit 
cell. Two types of symmetry of the molecules' positions 
inside such a plane lattice are the subject of our study, 
notably 

(a) Monoclinic or triclinic symmetry. This model mim- 
ics the (a, 6)-plane of polyacenes. We use some crystal- 
lographic and spectroscopic data for the crystals of an- 
thracene, tetracene, and naphthalene in our 2D models. 
The polyacene crystals exhibit layered structure which is 
better pronounced in the crystals of increasing number 
of the benzene rings. Our hypothesis of simple lattice 
neglects the effect of Davydov splitting (that effect has 
been treated, e.g., in the paper by Warns et al}'^ and in 
Lalov et al}^). The model of simple lattice allows us to 



obtain analytical and relatively simple results. 

(b) Hexagonal symmetry. The 2D sheets of hexagonal 
symmetry are actual nowadays not only because of their 
connection to the study of graphcnc but also in the crys- 
tal engineering of layered hexagonal structures (see, e.g., 
Thalladi et al}^). 

In the present paper we follow the vibronic ap- 
proach developed and successfully applied in our previous 
papers, ^"^^^^ in which one-dimensional models have been 
considered. Our present 2D study seems to be more real- 
istic, especially for the case of polyacenes. The 2D mod- 
els, however, limit the opportunity to treat the mixing of 
FE and CTEs. Whereas in ID models this mixing can be 
described using three exciton branches (one of FE plus 
two of CTEs) , the 2D models must include more than one 
decade of branches (see Petclenz et al.,^^ where the num- 
ber of mixed exciton branches is 14). The intention to 
simulate the details of the vibronic spectra and to inter- 
pret their structure is the reason to limit ourselves with 
the model of vibronic spectra of FE only in the simple lat- 
tice with one molecule per unit cell. Moreover, according 
to nowadays concept the lowest singlet exciton in poly- 
acenes is the Frenkel exciton and this is an argument for 
the treatment of FEs and their vibronics without mixing 
with CTEs. 

In calculating the linear optical susceptibility, Xj we 
use the formalism of the Green functions at T and 
the nearest neighbor approximation. In this way, we ex- 
press X s-iid calculate the linear absorption in terms of 
the complete elliptic integrals of the first kind 
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The outline of the paper is the following: in the sec- 
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ond section wc introduce the Hamiltonian and calculate 
the hnear optical susceptibihty in the range of one- and 
two-phonon vibronic spectra of a simple 2D lattice. In 
Sec. 3 those general expressions are specified for mono- 
clinic and triclinic 2D lattices and using crystallographic 
and spectroscopic data for polyacenes we make simula- 
tions of their excitonic and vibronic spectra. Our calcula- 
tions concern the vibronics of anthracene, tetracene, and 
naphthalene with intramolecular vibration at frequency 
of 1400 cm~^, but for the last crystal we simulate also 
the well studied vibronics with vibrations at frequency of 
702 cm~^. Section 4 deals with the vibronics in a hexag- 
onal 2D lattice with one molecule per unit cell. Two 
types of FEs are the subject of our studies of the exci- 
tonic density of states (DOS) and of linear absorption, 
namely (i) non-degenerate FEs whose transition dipole 
moment is perpendicular to the plane of the lattice, and 
(ii) degenerate FEs whose transition dipole moments are 
parallel to the layer. Section 5 summarizes our findings 
and contains some conclusions. In the Appendix, the 
Hamiltonian of degenerate FEs in the case of a hexago- 
nal lattice is established to split into two fully identical 
Hamiltonians of left and right FEs. 

II. HAMILTONIAN AND LINEAR OPTICAL 
SUSCEPTIBILITY IN A SIMPLE 2D LATTICE 

We consider the excitonic and vibronic excitations in 
a 2D lattice (a, b) with one molecule per unit cell. In 
the nearest neighbor approximation the FE part of the 
Hamiltonian reads 

Hf=J2 ^FB + Br, + Wnrr.B+Bm, (2) 

where i?„ (B+ ) is the operator of annihilation (creation) 
of the electronic excitation on the molecule n, Ep is the 
excitation energy of a molecule in the layers, and Wnm is 
the transfer integral of FE between molecules n and m. 
One mode of the intramolecular vibration at frequency 
Wo is supposed to be coupled with the FE and a„ is the 
annihilation operator of one vibrational quantum on the 
molecule n. Then the phonon part of the Hamiltonian 
can be written down as 

Hph = ^ hhJno^an- (3) 

n 

We suppose both linear and quadratic FE-phonon 
coupling^^^ with a Hamiltonian in the form 

^cx-ph = ^ ^fi^oB^ Bn {an + a^) 

n 

-h ^nAa;S+ B„a+a„, (4) 

n 

where ^ is a dimensionless parameter characterizing the 
linear FE-phonon coupling and Ao; is the change of the 



vibrational frequency of a molecule with electronic exci- 
tation on it (quadratic coupling). The three parts (2), 
(3), and (4) of the Hamiltonian can be transformed with 
the goal of eliminating the linear coupling by using the 
canonical transformation^'^ 

Hi=eMQ)Hexp{~Q), (5) 

where 

Q = Y.^^nBn{a+~a„). (6) 

In the transformed Hamiltonian the linear coupling 
terms are absent but the operators Bn are replaced by 

Vn = exp{Q)Bn exp(-g) (7) 

(for more details see Lalov and Zhelyazkov^^). In the 
momentum space (fca, kb) the operator (2) is transformed 
into 

Hf=Y.[Ep + W {ka,h)] V+^,Vk^,k,, (8) 

ha Mb 

where the matrix element W {ka, fcf,) of the intermolecu- 
lar transfer depends on the symmetry and on the strength 
of the intcrmolccular interaction. 

Wc use the following formulas^'"'^^ in calculating the 
linear optical susceptibility 

X^J = Inn I - ^ [^^J (^ + «e) + i-uj + ie)]^ (9) 
with 

$,,(i) = -ieit){o\p,it)p,io) + P,{t)P,{o)\o), (10) 

where V is the volume of the crystal (layer) , and P is the 
operator of the transition dipole moment. The Green 
functions (10) have been calculated as an average over 
only the ground state |0) taking into account the large 
values of Ep and that Hujo ^ ksT. In the expression of 
the operator P wc preserve the transition moment of the 
FE only that reads 

Pf=5]Pf(K + 1^+). (U) 

n 

Here Pp is the transition dipole moment of the electronic 
excitation on the molecule n. The calculations of % are 
reduced to calculating the following Green functions: 

G(°)(t) = -z0(i)(O|K(OV^+(O)|G). (12) 

In calculating their Fourier transforms in momentum 
space, G'^°'(0, 0), wc obtain the following equation: 

[huj- Ef ~W{0, 0)] G(°) (0, 0) = 1 + huJaS'^^^ , (13) 
where 

Wa=?Wi, UJi^UjQ + Auj, (14) 
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and 5^^^ is the sum on the whole 2D Brillouin zone 
(fed, fcf,) of the Fourier transforms of following Green func- 
tions: 

Gi%it) = -~te{t){OMt)Vn{t)V+mO). (15) 

In calculating the Green functions (15) we ob- 
tain functions G'^^ with two phonon operators 



{0\an{t)ajn{t)Vn{t) ■ ■ ■ and correspondingly one ladder 
with more complicated Green functions appears. Namely 
calculating them we use the vibronic approach^° in which 
the transfer terms W {ka, k^) are taken in one step of the 
ladder only. 

If we are interested in one-phonon vibronic spectra, the 
following equation of the component G^^^ {ka, kb) is valid 



{h [oj-cuo- ^^r(l)] -Ef-W ika,kb)} G(i) (fc„ kb) = hwaG^''^ (0, 0) + a^(i), (16) 

where 

a ^ h[no^r{2) - nr{l) + Alj] , (17) 
and llo,r(2) and 17^(1) ai'c the following continuous fractions: 

f^o,r(2) = ^-^ , (18) 

uj — Ep /ft — 2a;i — 



uj — E-p jh — 3aji 



UJ — Ep /Ti — 2uj{) — Aw — 



UJ — E-p /h — Sujq — 2Auj 
After some algebra we find the following expression of the Fourier transform G(°)(0, 0) 

G(o)(0,0) = — 



Huj - Ep ^WiO,0) 



1-aTi 



where 



Ti = ^ {ft [W - wo - f^r(l)] -Ep-W {ka, kb)}-' . 



ka -kh 



Finally we find the following expression of the linear optical susceptibility: 

X.. = -A:^g(")(0,0), 



(19) 



(20) 



(21) 



(22) 



in which A depends upon the units, v is the volume occupied by one molecule, and the a;-axis is directed along the 
direction of the transition dipole moment Pp. In the region of the two-phonon vibronic spectra the same procedure 
yields the expression^" 
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huj-Ep- W{Q, 0) 



huj-i 



(23) 



UJ — Ep/h — uji — 



2ujIT2 
1 — a2T2 



where 



02 = 2hAuj + hujl 



— Ep /ft — 3wi — 



4uji 



■- 0,r(2) 



(24) 



Ep/h — 4wi 
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a.(2) = 2 , (25) 

LO — E-p /h — SuJq — Aid 



3^2 

Lj - Ep/h - iujQ - 2Auj 



and 

^ {h[uj - 2u;o - nr{2)] - Ef - W {ka,h)y\ (26) 

III. TWO-DIMENSIONAL LATTICE OF MONOCLINIC AND TRICLINIC SYMMETRY 

In this case, the molecules create a plane rhombohedral or rectangular network with two nearest neighbors at 
a-direction and other two nearest neighbors at 6-direction. We denote by Wa/2 the transfer integral between the 
neighbors in the a-direction and by Wb/2 the transfer integral in the 6-direction. The quantity W (kajkh) can read as 

W {ka, h) = Wa cos [kaO) + Wb COS (fefefe) , (27) 

and we find the following expressions of the quantity Ti (and T2 as well) : 

sin B 

Ti = , ^ J, (28 

where /? is the angle between the a* and b* axes in the reciprocal space. To simplify notation we introduce also the 
following quantities: 

Pi^uj - Ep/h- ujQ-Vlr{l) and ^ UJ - Ep / h - 2ujq - D.r{2) , (29) 



S/3; 



— for onc-phonon vibronic. 



\Wb\ 

for two-phonon vibronic, 

\Wb\ 



(30) 



Wa 



ki-d^^J^^ k2-d . (31) 



Wb 

Then the integral / in Eq. (28) can be calculated (assuming \ Wb\ > \Wa\) as follows: 

/ = -kiK{ki) for Re(i) < -1 -p, (32a) 
I = -K{l/ki) ~iK{k2) for - 1 < Re(t) < -1 +p, (32b) 

/ = -^K{l/k2) for -l+p< Rc(t) < 1 - J5, (32c) 

I = K{l/ki) -iK{k2) for 1 -p < Re(t) < 1 +39, (32d) 
/ = kiK{ki) for Re(t) > 1 + p. (32e) 

The same expressions are valid for T2 but with substituting /32 for /3i . 

We note here that the quantity Ti expresses the excitonic density of states (DOS) at wo = 0, ^ = 0, and Aw = 0. 
Then we need the Fourier transform G^*^^ (fcj,, kb) of the retarding Green function (12) which enters the expression of 
DOS, p{lo), in the excitonic band: 

p{uj) = -Trim ^ G(°' (fc, w) at = 0, ^ = 0, Aw = 0, u; ^ u; + i(5 for (5 ^ 0. (33) 

k 

Finally, it is easy to obtain the following relationship: 

p{uj) = -TTlmTi (wo = 0, C = 0, Aw = 0) when <5 ^ 0. (34) 
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A. Simulations of excitonic DOS and vibronic spectra of 2D models of polyacenes 

Naphthalene and anthracene crystals are of monoclinic symmetry and their crystallographic axes a and b are 
orthogonal as are a* and b* of their reciprocal lattices. Tetracene crystals are triclinic and the angle ao between 
the a and b axes is approximately equal to 100° but the a* and b* axes in the reciprocal lattice are almost perfectly 
orthogonal. 

In the spectra of the three crystals we are largely interested in the intensive and widely studied electronic transition 
whose transition electric dipole moment Pf is directed along the M-axis of the molecule I (see Campbell et al.^^). 
In our model we calculate the transition integrals Wa and Wb due to the nearest-neighbor dipole-dipole interaction 
of the equal dipole moments Pp of each molecule of the 2D (a, 6)-lattice. The distance between neighbor dipoles is 
taken to coincide with the lattice parameters a and b. 

Table I contains some crystallographic data of polyacenes crystals^'^'^''^ which have been used in the simulations. 
Table II contains spectroscopic data necessary for our simulations. We take those data from Schlosser and Philpott^^'^^ 



TABLE I. Crystallographic data of polyacenes. 



Crystal 


Symmetry 


a (A) 


6(A) 


a = [a, b) 






Naphthalene 


Monoclinic 


8.24 


6.00 


90° 


29.5° 


^69.6° 


Anthracene 


Monoclinic 


8.56 


6.04 


90° 


26.6° 


71.3° 


Tetracene 


Triclinic 


7.90 


6.03 


100.3° 


30.1° 


69.2° 


" Angle is between the M and b axes, 


while the angle t/jM 


is between the M and 


ai ± b-axis (ai is lying 


in the [a, ii)-plane 


)■ 


to calculate the transfer integrals Wa 


and Wb as the nearest neighbor dipole-dipole interaction. Recall that in 


Table II 


fi = Pp/e is the dipole transition moment measured in angstroms. In naphthalene the lowest excitonic transition is very 




TABLE II. Spectroscopic data for excitonic and vibronic spectra of polyacenes 


24 




Crystal 


Ef (oV) 


fiojo (eV) 


M = Pp/e (A) 


Wa (eV) 


Wb (eV) 


e 


Naphthalene I 


4.33 


0.1772 


0.54 


0.0097 


-0.049 


1.67 


Naphthalene II 


3.87 


0.0942 


0.104 


-1 X 10"^ 


-0.011 


0.4 


Anthracene 


3.11 


0.1735 


0.61 


0.012 


-0.068 


0.9753 


Tetracene 


2.615 


0.1772 


0.69 


0.012 


-0.078 


1.207 



weak (transition II in Table II) with an oscillator strength of the polarization parallel to the 6-axis, fb « 4x 10-^ and 
respectively for polarization parallel to the a-axis, /a ~ 2 x 10~^ (sec Chap. 2 in Broude et al.^). The calculated values 
of Wa and Wb are correspondingly much lower in comparison to other transfer integrals. The previous spectroscopic 
studies of molecular and crystal spectra exposed a significant shift, hAui = —0.007 eV, of the vibrational frequency 
in an excited molecule. 



200 - 



150 - 




3 3.05 3.1 3.15 3.2 



Fico (eV) 

FIG. 1. (Color online) Excitonic DOS in the case of dispersion (27) in the excitonic band. Ef — 3.1 eV, Wa = 0.012 eV, 
Wb = —0.068 eV. The red curve is calculated at hS = 0, and the green one at M = 1 x 10~^ eV. 
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In the following we use relative units and suppose that the product APp/v = 1. We add an imaginary part, iS, to 
uj which expresses the excitonic damping and perform the calculations with hS = 1 x 10"'^ eV. The linear absorption 
coefficient is calculated as imaginary part of the component Xbb (for electromagnetic waves of polarization parallel to 
the 6-axis) 

X6fc = -cos2 7^MG(°)(0,0). (35) 

The excitonic DOS calculated by using formula (33) is shown in Fig. 1. The well-known Van Hove singularities^^ 
for the case of 2D models are exhibited. The function of DOS has non-zero values at the edges of the excitonic band 
and it manifests two singular points of logarithmic behavior. The final excitonic damping (the green curve) preserves 
the course of DOS but it makes the singularities softer. 




FIG. 2. (Color online) Vibronic spectra of the anthracene 2D model. The excitonic peak is fitted at Ef ~ 3.11 eV, Wa = 
0.012 eV, Wb = -0.068 eV. The red curve is calculated at = 0.6, the green at = 0.8, the blue at = 0.9753, and 
the black at — 1.2, all curves for Acj — and fkoo = 0.1735 eV. Clip (a): general picture of the excitonic peak and the 
absorption curves of one-phonon vibronics (formula (20)). Clip (b): absorption curves near Ep + hujo (the first replica). Clip 
(c): absorption curves near Ef + 2hujo (second vibronics, formula (23)). 



70 




3.43 3.44 3.45 3.46 3.47 3.48 3.49 3.5 



ho) (eV) 

FIG. 3. (Color online) Two-phonon vibronic spectra of the anthracene 2D model. The red curve has been calculated for 
^2 = 0.8, hAuj = 0; green curve for ^2 = 0.8, ftAoj = -0.01 eV; blue curve for = 0.9735, hAuj = 0; and the black curve for 
= 0.9735, hAuj = -0.01 eV. 

The linear absorption simulations of the anthracene 2D model calculated on using formulas (20) and (23) yield the 
curves in Fig. 2. In our calculations, the value of Ep is a free parameter (see Ref. 20) and it is fitted to obtain the 
excitonic peak at 3.11 eV. For small values of the linear exciton-phonon coupling (red and green curves) the absorption 
curves in Fig. 2(b) are wide and correspond to the many-particle (MP) states. The half-width of the black curve is 
approximately of 1 x 10""^ eV and it is Lorcntzian which corresponds to the boimd (one-particle) exciton-phonon 
state. The blue curve is calculated for values of ^ close to the widely used values of anthracene and it corresponds to 
a quasi-bound state inside the MP band near its minimum. The jumps in the absorption curves are associated with 
the singular points of DOS. 

All curves in Fig. 2(c), calculated for the second vibronics Ep + 2haja, exhibit behavior of MP states. The position 
of the MP bands can be found using the same formulas (20) and (23) but at 5 = 0. In this case, one-particle states 
do not manifest themselves in the linear absorption spectra but MP bands appear due to the imaginary parts in 



7 




4.3 4.4 4.5 4.6 4.7 4.8 4.9 4.48 4,485 4.49 4.495 4.5 4.505 4.51 4.515 4.52 4.525 4.53 4.66 4.665 4,67 4.675 4.68 4.685 4.69 4.695 4.7 4,705 4.71 



lito (eV) hco (eV) Tim (eV) 



FIG. 4. Linear absorption spectra of the naphthalene I model. Ep = 4.33 eV; tiojo = 0.1772 eV; = 1.67; Acj = 0; 
Wa = 0.0097 eV, and Wb ~ —0.049 eV. Clip (a): general picture of the vibronic series calculated on using formula (20). Clip 
(b): absorption curves near Ef + huo. Clip (c): absorption curves near Ef + 2&Jo, calculated by means of formula (23). 




FIG. 5. (Color online) Linear absorption of the lowest excitonic peak and its vibronic replica of naphthalene II (see Table II). 
Ef = 3.87 eV; hujo = 0.0942 eV; = 0.4; hAuj = -0.0072 eV. The green curve has been calculated for Wb = -0.011 eV, 
Wa = -1 X lO"'^ eV; for red curve Wb = 0.011 eV, = 1 x 10"^ eV. Clip (a): general picture. Clip (b): the region of the 
first vibronic spectra Ef + hujQ. A factor of 0.04 decreases the values of the absorption in comparison to Fig. 4. 




2.6 2.65 2.7 2.75 2.8 2.85 2.9 2.95 




rico (eV) 



fico (eV) 



FIG. 6. (Color online) Linear absorption of the tetracene 2D model. Ef = 2.615 eV; hujo = 0.1772 eV; = 1.207; Wa = 
0.012 eV, Wb = —0.0781 eV. The red curve corresponds to Au = 0, and the green curve to hAuj = —0.02 eV. Clip (a): 
absorption curves calculated by using formula (20). Clip (b); absorption curves near Ef + 2/kJo calculated with the help of 
formula (23). 



formulas (32). The calculations a.t 6 — show that even the most narrow band for = 0.6 is positioned inside the 
MP band. It is curious to note the inverse situation in vibronic spectra with two phonons compared to Fig. 2(b). The 
increasing linear excitation makes the absorption curves wider (in one-phonon vibronic spectra it binds the FE and 
phonon), however, shifts them in the region of lower frequencies like curves in Fig. 2(b). The impact of the quadratic 
exciton-phonon coupling expressed by Aw is more pronounced in the two-phonon vibronic region (see Fig. 3). The 
shift of the absorption curves to lower frequencies is even more than 2hAuj (compare red and green curves, blue 
and black curves) but the quadratic coupling transfers intensities from two-phonon vibronic to one-phonon vibronic 
spectra. 
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The vibronic series of the naphthalene I model (see Table II) is shown in Fig. 4. Obviously the first vibronic 
replica (look at Fig. 4(b)) is dominated by the strong one-particle maximum accompanied by a relatively weak MP 
band in the region 4.505-4.53 eV. The situation near the second vibronic (see Fig. 4(c)) is very similar, however, the 
one-particle maximum and the MP band have comparable absorption intensities. 

The excitonic maximum, E-p = 3.87 eV, and its vibronic replica have been intensively investigated^'^'^ despite 
of their weak intensities. Our simulations shown in Fig. 5 confirm the results obtained by using an ID model (see 
Ref. 20). In the case of negative values of the transfer integrals Wa and Wb the absorption spectrum consists of 
onc-particle maximum and a very weak MP band (green curve in Fig. 5). In the case of positive transfer integrals 
(red curve) the linear absorption demonstrates the transfer of unbound FE and phonon. 

The linear absorption of the tetracene 2D model presented in Fig. 6 exhibits one-particle maximum near + Hujq 
and a MP band near the second vibronic replica. In polyacenes, the vibration of quantum hioo « 0.17 eV is an example 
of presumably linear FE-phonon coupling and thus the red curves in Fig. 6 are more realistic. The hypothetical strong 
quadratic coupling causes the binding of FE and phonons and one-particle states only occur in the linear absorption. 



IV. TWO-DIMENSIONAL LATTICE OF HEXAGONAL SYMMETRY 



Each molecule in a plane hexagonal lattice with one molecule per unit cell is surrounded by six neighbors spaced 
at a distance a. We us the unit cell with vectors a and b of displacement from molecule to molecules 1 and 2 and 
the angle ao between a and b is 60° (|a| = |b| = a). The angle /3 between the a* and b* axes of reciprocal space is 
equal to 120°. 

Two types of FEs will be considered in the following text: (i) non-degenerate FEs of transition dipole moment 
perpendicular to the plane of the hexagonal layer, and (ii) degenerate FEs whose transition dipole moments are 
parallel to the layer. 

The transfer integral Vi between the molecule and its six nearest neighbors can be calculated using the formula 
for the potential energy of two molecular transition dipole moments pi and p2 whose centers are spaced by a vector r 

W = (Pi-P2)^^-3(Pi-r)(piT) 

where Eq = 8.8542 x 10"-^^ Fm~-^ is the electric constant. 

In the case of non-degenerate FEs the transition dipole moments |pi| = |p2| = p are perpendicular to the vector r 
and we obtain 

In the case of degenerate FEs their transition dipole moments directed along the perpendicular axes x and y inside 
the plane of the layers are equal: hx = hy = h (see the Appendix). Using the same formula (36) one obtains the 
transfer integrals of both left and right FEs 

Vi = -^-^ < 0. (38) 

Irrespective of the different signs of the transfer integrals Vi of non-degenerate and degenerate FEs one obtains 

W (ka, h) = 2Vi {cos(7fca) + cos(7fc&) -I- cos [7 (fca + h)]} , (39) 
where 7 = aV3/2. We introduce the quantities 

s=^, d=V3~2-s, (40) 



Then the sums Ti and T2 (see formulas (21) and (26)) can be expressed using the integral /: 
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whose new values now are given by 



2^/d 
1 



/ = - 



7 = 



2Vd 
k2 



2Vd 



A'(fci) for Re(s) < -3, 

[K{l/ki) -isgn{Vi)K{l/k2)] for - 3 < Re(s) < 1, 



3 - s2 - 2rf 



3 - s2 + 2fi 



iSgn(V^i)K(fc2) 



1 



1 



for 1 < Re(s) < 3/2, 

for Re(s) > 3/2. 



(43a) 
(43b) 

(43c) 
(43d) 



The excitonic DOS in the case of excitonic dispersion (39) can be calculated with the help of formula (34). In 
the simulations of DOS and linear optical absorption we put the excitonic vibrational parameters typical for organic 
solids, more specifically 

hujo ^ 0.17 eV, 
= 0.64, or 1, or 1.44, 



hAuj = 



- 0.01 eV. 



In estimating the transfer integral Vi we suppose a dipolc-dipole approximation between two transition dipoles 
p (or h/^/2) = 3 D situated at a distance a — 8 x lO^^*' m. Then the calculations of their potential energy give 
(approximately) 

(a) Vi = 0.011 eV for non-degenerate FEs. We fit the excitonic level at Ep — 3.27 eV. 

(b) Vi = -0.011 eV for degenerate FEs with Ep = 3.06 eV. 

The excitonic DOS in the case of dispersion low (39) is presented in Fig. 7. Only one saddle point that corresponds 
to Re(s) = 1 (see formulas (40)~(43)) exists in the excitonic band. In the center of the Brillouin zone, ka = kb = 0, 
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FIG. 7. (Color online) Excitonic DOS in the case of dispersion described by formula (39). For non-degenerate FEs Ep = 3.27 eV, 
Vi = 0.011 eV, while for degenerate FEs Ef = 3.06 eV and Vi = -0.011 eV. 




the values of DOS at Re(s) = —3 are relatively low. 



A. The case of degenerate Frenkel excitons (Vi < 0) 

In Fig. 8(a) one sees the calculated absorption curves near the excitonic peak and one-phonon vibronic spectra. In 
Fig. 8(b) are shown the lowest parts of the same dispersion curves calculated at S = (no excitonic damping) in the 
bands of the many-particle FE-phonon states. We emphasize again that the absorption at S — manifests itself in 
the MP bands only. The points of a sharp change in the three absorption curves correspond to the singular point 
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FIG. 8. (Color online) Linear absorption in a simple hexagonal 2D lattice for Ef — 3.06 eV, hAuj = 0.17 eV, hAuj — 0, 
Vi = —0.011 eV. The green curve has been calculated at ^ = 0.8, the red for ^ = 1, and the blue for ^ = 1.2. Clip (a): excitonic 
peak ant its vibronic replicas at hS — 1 x 10"'' eV. Clip (b): Details of the absorption curves in the MP bands at 5 = 0. 



Rc(s) = 1. Both panels of Fig. 8 show that the green curve for f = 0.8 describes absorption in the MP band, the 
red curve = 1) corresponds to a quasi-one particle state above the minimum of the MP band while the blue curve 
~ 1.2) describes one-particle (bound) FE-phonon state below the MP band. 

The impact of the quadratic FE-phonon coupling {hS ^ 0) on the absorption curves is clearly seen in Fig. 9. The 
absorption maxima are shifted approximately at HAijj below their positions in Fig. 8(a) and the maximum of the red 

700 1—1 — , 1 1 1 1 1 

600 - 



500 - 




3.05 3.1 3.15 3.2 3.25 



hco (eV) 

FIG. 9. (Color online) Linear absorption as in Fig. 8 but for hAcj — —0.01 eV. 

curve (^ = 1) describes also a one-particle state. 

Figure 10 illustrates the absorption near the second vibronic replica Ep + 2hujQ. The linear FE-phonon coupling 
is not sufficient to bind the FE and two phonons and the three absorption curves in Fig. 10(a) describe MP states. 
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FIG. 10. (Color online) Linear absorption as in Fig. 8 near the second replica, E-p -f 2kujo. Clip (a): the case hAui = 0. Clip 
(b): the case hAuj = -0.01 eV. 
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The calculations performed at 5 = 0, similar to those in Fig. 8(b), confirm the position of the maximum on the green 
curve inside the MP band irrespective of its narrow width. The quadratic coupling shifts the two-phonon vibronic 
maxima in the direction of lower energy — see Fig. 10(b). Unlike the case of linear coupling (look at Fig. 10(a)) the 
maximum of the blue curve (^ = 1.2) corresponds to bound states while the maximum of the red curve (^ = 1) is 
associated with a quasi-bound state. 

Our studies of the linear absorption show that both coupling mechanisms decrease the frequencies of the absorption 
maxima of the degenerate FE-phonon states (Vi < 0). 



B. The case of non-degenerate Frenkel excitons (Vi > 0) 



The picture of the absorption curves at Vi > seems to be a little surprising (see Fig. 11). The absorption maxima 
appear near the maximum of the MP band (green curve) and even above the PM bands (red and blue curves). The 
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FIG. 11. (Color online) Linear absorption in a 2D model of hexagonal lattice — the case of non-degenerate EEs — excitonic peak 
and first vibronic spectra. Here Ef = 3.27 eV, hujo = 0.17 eV, hS = 1 x 10~^ eV, Vi = 0.011 eV. The green curves have been 
calculated at ^ = 0.8, the red curves for ^ = 1, and blue ones for ^ = 1.2. Clip (a): the case Aw = 0. Clip (b); the case 
hAu = -0.01 eV. 



distance between the excitonic peak (E-p — 3.27 eV) and onc-phonon vibronic maxima is bigger than hwQ = 0.17 eV. 
Hence, the linear FE-phonon coupling at > repulses the vibronic levels and increases their energy. In contrast, the 
quadratic coupling (see Fig. 11(b)) decreases the values of the vibronic levels and one-particle bound states above the 
MP bands in Fig. 11(a) are transformed into unbound states inside the MP bands (red and blue curves in Fig. 11(b)). 
The three absorption curves in Fig. 12 correspond to the MP bands. The absorption at ^ = 1.2 (blue curve) 
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EIG. 12. (Color online) Linear absorption spectra as in Eig. 11 near the second vibronic replica, Ef + 2hwo, calculated by 
means of formula (23). Clip (a): the case fiAo; = 0. Clip (b): the case h/S.uj = —0.01 eV. 



exhibits also the existence of one-particle state — see the high maximum above the energy of 3.705 eV in Fig. 12(a). 
Like the one-phonon vibronic spectra (see Fig. 11) the simultaneous action of the linear and the quadratic couplings 
(look at Fig. 12(b)) shifts in opposite directions the two-phonon vibronic spectra. We stress on the difference between 
absorption curves in Fig. 11(b) (calculated by using expression (20)) and the absorption in the same region 3.45- 
3.50 cV in Fig. 12(b) (calculated with the help of formula (23)). Formula (20) concerns onc-phonon vibronic spectra 
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(see Fig. 11) and is not valid for the one-phonon maxima calculated using the expression (23) that describes the 
two-phonon vibronic spectra. 

The most interesting result of our 2D model in the case of a FE polarized perpendicularly to the layer is the opposite 
action of the two coupling mechanisms on the position of the vibronic maxima. 

I 



V. CONCLUSION 

In the present paper, we investigate the excitonic and 
vibronic spectra of FEs manifesting themselves in 2D 
plane lattices with one molecule per unit cell of the fol- 
lowing types: (i) monoclinic or triclinic lattice which 
mimics the (a, b) plane of polyacenes, and (ii) simple 
hexagonal lattice. We study the exciton dispersion and 
the excitonic DOS in the nearest neighbor approximation 
as well as calculate the linear optical susceptibility and 
the linear absorption spectra in the range of a FE and 
its vibronics with one and two quanta of intramolecular 
vibration coupled with the FE by linear and quadratic 
couplings. The model of a lattice with one molecule per 
unit cell yields relatively simple expressions with com- 
plete elliptic integrals of the first kind. Our approach 
allows us to find the positions of the MP (unbound) FE~ 
phonon bands and of the onc-particlc (bound) states out- 
side those bands. The Van Hove singularities appear in 
the DOS of the 2D models and they affect the linear ab- 
sorption spectra. 

In the case of a hexagonal 2D lattice our studies con- 
cern two types of dipole-active FEs: (i) non-degenerate 
FEs of transition dipole moment perpendicular to the 
layer and of positive transfer integral, W; (ii) degenerate 
FEs of transition dipole moments parallel to the layer 
and of negative transfer integrals. In the last case, we 
establish a splitting of the Hamiltonian into two fully 
identical Hamiltonians which do not mix in the dipole 
approximation and describe two types of FEs of different 
hiralitics (left and right, respectively). The simulations of 
the vibronic spectra exhibit a opposite impact of the lin- 
ear FE-phonon coupling on the positions of the vibronic 
maxima: in case (i) of non-degenerate FEs it repulses 
those maxima to the upper boundary of the MP bands 
or above these bands. In case (ii) of degenerate FEs the 
linear coupling decreases the energy of the vibronic max- 
ima. Since the quadratic coupling in the usual case of 
HAlo < decreases this energy, the simultaneous action 
of both coupling mechanisms can have a cumulative effect 
in case (ii) or a compensating one in case (i). 

The 2D simulations of the vibronic spectra of naphtha- 
lene and anthracene in the present paper agree well with 
the ID simulation in Ref. 20, however, new peculiarities 
in the absorption near the singular points could appear. 
Especially inside the two-phonon MP bands narrow vi- 
bronic maxima can be observed even at small values of 
the linear coupling constant ^ (see Figs. 2(c), 10(a), and 
12(b)). 

The intriguing models of layered polyacenes and of 
graphene, both with two molecules per unit cell can not 
be described using the relatively simple expressions of the 



present paper — one needs the usage of elliptic integrals of 
the third kind for an adequate modeling. Nevertheless, 
many conclusions for the absorption inside and outside 
the MP bands are still valid. The studies of the excitonic 
dispersion and the structure of the vibronic spectra can 
be useful in the modeling of the excitonic transfer and 
the exciton-phonon coupling manifesting themselves in 
other linear and nonlinear phenomena. 

Appendix: Hamiltonian of degenerate Frenkel 
excitons 

The dipole-active FEs whose transition dipole mo- 
ments are parallel to the hexagonal layer are two-fold 
degenerate. The same two-fold degenerate electronic ex- 
citations must exist in each molecule of the monolayer. 
In the lack of rotational symmetry of molecules in the 
(x, y)-plane of the monolayer it is impossible to ensure 
optical isotropy inside that plane. 

We choose the axes x and y to be perpendicular to 
each other and parallel to the monolayer. The Hamilto- 
nian and the transition dipole moment of each molecule 
contain the following components: 

i7„,oi = Ef {B+B,. + B+By) , (A.l) 

Pmoi = h [{B+ + B,) i + {B+ + By) y] , (A.2) 

where B^ and By arc the operators of annihilation of the 
electronic excitations in a molecule with transition dipole 
moment directed along the x/y-axis, and x and y are the 
corresponding unit vectors. We introduce the following 
operators: 

K = {Bt + ^B+) /V2, Bi = (B, - iBy) /^/2, 
B+ = {B+ - iB+) /V2, Br = {B, + iBy) /V2. 

The usual boson commutation rules are valid for the 
operators (A. 3). Thus, operators (A.l) and (A.2) take 
the forms 

^mol = Ef {B+Bi + B+Br) , (A.4) 

Anoi = {h/V2) [{x - ly) {B+ + Br) 

+ {x + ly) {B+ + Bi)] . (A.5) 

The point group of symmetry of the hexagonal plane 
monolayer can be Ce, C3/1, or D^h depending on the sym- 
metry of the molecules. In all cases the two-dimensional 
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{x, ?/)-representation of the point group can be pre- 
sented by the foUowing two one-dimensional represen- 
tations whose direct product is a totally symmetrical 
unit-representation: (i) left (l) with molecular transition 
dipole moment {h/V2){x — iij), and (ii) right (r) with 
transition dipole moment {h/^/2){x + iy). 

The non-vanishing terms in the crystal Hamiltonian 
must be invariant for all operations of symmetry of the 
crystal and the binary terms of the energy operator and 
that of the intcrmolccular interaction must be a product 
of the two different representations, r and I, respectively. 
In the Heitler-London approximation (sec Davydov^ and 
Agranovch-'^^) we neglect the terms 



and the only non- vanishing transfer terms are of the type 



In dipole-dipole approximation the corresponding trans- 
fer integral Vi can be calculated using formula (39) for 
the two dipoles 

Pi = {h/V2){x - iy) and p2 = {h/V2){x + iy). 

In this way, one obtains formula (38). The transfer in- 
tegrals, Vi, for the excitations I and r are equal and 



the excitonic Hamiltonian, as well as operator (4) of 
the exciton-phonon coupling, split into two independent 
parts / and r which do not mix, notably 



^""'^(t,r);m^(','-);" 



(A.6) 



Hence, the Hamiltonian Hi can be divided into two 
independent Hamiltonians Hi^i and Hi r fully analogous 
to Eq. (8). In calculating the linear optical susceptibility, 
X, we apply the approach used in Sec. II for the transition 
dipole moment (A. 5). Since FEs of type r and I do not 
mix, the non- vanishing Green functions of the type (12) 
in the Hamiltonian must have equal chirality r or Z of the 
operators V{t) and V^{0). The final expression for the 
components 



Xyy 



(A.7) 



is fully identical to (22) but with h instead of Pp and new 
values of the transfer integrals Vi which are negative — 
compare expressions (37) and (38). 



* izh@phys.uni-sofia.bg 

^ W. T. Simpson and D. L. Peterson, J. Chem. Phys. 26, 
588 (1957). 

^ A. S. Davydov, Theory of Molecular Excitons (Plenum 

Press, New York, 1971). 
3 R. E. Merrifield, J. Chem. Phys. 40, 445 (1964). 
" E. I. Rashba, Sov. Pliys. JETF 23, 708 (1966). 
^ M. R. Philpott, J. Cliem. Pliys. 47, 2534 (1967). 
^ E. I. Rashba, Sov. Phys. JETF 27, 292 (1968). 

M. R. Philpott, J. Chem. Phys. 55, 2039 (1971). 

* E. F. Sheka, Sov. Phys. Usp. 14, 484 (1972). 

® V. L. Broude, E. I. Rashba, and E. F. Sheka, Spectroscopy 
of Molecular Excitons (Springer, Berlin, 1985). 
V. M. Agranovich, Excitations in Organic Solids (Oxford 
University Press, New York, 2009). 
" P. Petelenz, M. Slawik, K. Yokoi, and M. Z. Zgierski, J. 
Chem. Phys. 105, 4427 (1996). 

M. H. Henessy, Z. G. Soos, R. A. Pascal Jr., and A. 
Girlando, Chem. Phys. 245,199 (1999). 
M. Hoffmann, K. Schmidt, T. Fritz, T. Hasche, V. M. 
Agranovich, and K. Leo, Chem. Phys. 258, 73 (2000). 
" M. Hoffmann and Z. G. Soos, Pliys. Rev. B 66, 024305 
(2002). 

I. J. Lalov, C. Supritz, and P. Reineker, Chem. Phys. 309, 



189 (2005). 

I. J. Lalov and I. Zhelyazkov, Phys. Rev. B 74, 035403 

(2006) . 

C. Warns, I. Lalov, and P. Reineker, Phys. Procedia 13, 
33 (2011), Selected Papers from 17th International Con- 
ference on Dynamical Processes in Excited States of Solids 
(DPCTO). 

I. J. Lalov, T. Hartmann, and P. Reineker — submitted. 
V. R. Thalladi, K. Panneerselvam, C. J. Carrell, H. L. Car- 
rell, and G. R. Desiraju, J. Chem. Soc, Chem. Commun. 
issue 3, 341 (1995), DOI: 10.1039/C39950000341. 
I. J. Lalov and I. Zhelyazkov, Phys. Rev. B 75, 245435 

(2007) . 

I. J. Lalov, C. Supritz, and P. Reineker, Chem. Phys. 352, 
1 (2008). 

I. J. Lalov, C. Warns, and P. Reineker, New J. Phys. 10, 
085006 (2008). 

R. B. Campbell, J. M. Robertson, and J. Trotter, Acta 
Cryst. 15, 289 (1962). 

D. W. Schlosser and M. R. Philpott, Chem. Phys. 49, 181 
(1980). 

D. W. Schlosser and M. R. Philpott, J. Chem. Phys. 77, 
1969 (1982). 

L. Van Hove, Phys. Rev. 89, 1189 (1953). 



